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£h Lars Onsager announced in 1949 that he and Bruria Kaufman had proved a simple 

formula for the spontaneous magnetization of the square-lattice Ising model, but did not 
publish their derivation. It was three years later when C. N. Yang published a derivation 
in Physical Review. In 1971 Onsager gave some clues to his and Kaufman's method, and 
there are copies of their correspondence in 1950 now available on the Web and elsewhere. 
Here we review how the calculation appears to have developed, and add a copy of a draft 
fNJ paper, almost certainly by Onsager and Kaufman, that obtains the result. 
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cn 1 Introduction 

o 

Onsager calculated the free energy of the two-dimensional square-lattice Ising model in 
1944. [1] He did this by showing that the transfer matrix is a product of two matrices, 
• — which together generate (by successive commutations) a finite-dimensional Lie algebra 

(the "Onsa ger algebra"). In 1949 Bruria Kaufman gave a simpler derivation [2] of this 
result, using anti-commuting spinor (free-fermion) operators, i.e. a Clifford algebra. 

Onsager was the Josiah Willard Gibbs Professor of Theoretical Chemistry at Yale 
University. Kaufman had recently completed her PhD at Columbia University in New 
York, and was a research associate at the Institute for Advanced Study in Princeton. 

Later that year, Kaufman and Onsager[3] went on to calculate some of the two-spin 
correlations. Let i label the columns of the square lattice (oriented in the usual manner, 
with axes horizontal and vertical), and j label the rows, as in Figure [l] Let the spin at 
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site be cry, with values +1 and —1. Then the total energy is 

E = -J^CTijCTij+l - j'^CTijCTi+lj 



'■J 



and the partition function is 



z = Y,^ E/KT 



the sum being over all values of all the spins, k being Boltzmann's constant and T the 
temperature. Onsager defines H,H',H* by 



H = J/kT , H' = J'/kT , e~ 2H * =tanhH . 
The specific heat diverges logarithmically at a critical temperature T c , where 

sinh(2 J/kT c ) sinh(2 J'/kT c ) = 1 . 



'1.11 



(1,2) 



J 



Figure 1: The square lattice. 



The correlation between the two spins at sites (1, 1) and is 

(7 

For the isotropic case H' = H, Kaufman and Onsager [3] give in their equation 43 the 
formula for the correlation between two spins in the same rowQ 

(0-1,10-1,1+y) = cosh 2 H* Aj - sinh 2 H* A-j (1.2) 



1 We have negated their E r , which makes them the same as those in Appendix A, and corrected what appears 
to be a sign error. It is now the same as III. 43 of the draft paper below. 
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Here A,- and A_,- are Toeplitz determinants: 



A, 



Si 

So 



S2 S3 

Si S2 



s, 



A_ 



E_i S_2 S_3 

So S_i S_2 

Sj_2 • 



Si 



Si-j 



(1.3) 



where 



and 



1 

2^ 



2tt 



tan 



sinh2.£f sinw 



Setting 
this implies 



coth 2H' — cosh 2ff cos uj 



a W(w) 



1 - cothF' e~ 2H+iuJ ){l - tanhil' e" 2 ^+ ia; ) 
1 - cothif' e~ 2H ' iul )(l - tanhF' e - 2 ^-^) 



1/2 



(1.4) 
(1.5) 

(1.6) 



(Equations (1.4), (1.6) follow from (89) of [T] and are true for the general case when H' , 



are not necessarily equal.) 

Kaufman and Onsager also give the formula for the correlation between spins in adja- 
cent rows. In particular, from their equations 17 and 20, we obtain 



(0-1,10-2,2) = cosh 2 if* Si +sinh 2 iTS_i . 

The long-range order, or spontaneous magnetization, can be defined as 

1/2 



lim (anan+j) 



(1.7) 



(1.8) 



It is expected to be zero for T above the critical temperature, and positive below it, as in 
Figure [2] 

Kaufman and Onsager were obviously very close to calculating Mq: all they needed to 
do was to evaluate Aj, A_^ in the limit j — > 00. Here I shall present the evidence that they 
devoted their attention to doing so, and indeed succeeded. They used two methods: the 
first is discussed in section 2 and Appendix A, the second in sections 3, 4, 5. Many years 
ago (probably in the early or mid-1990's), John Stephenson (then in Edmonton, Canada) 
sent the author a photocopy of an eight-page typescript, bearing the names Onsager and 
Kaufman, that deals with the topic. Stephenson had copied it about 1965 in Adelaide, 



3 



Australia, from a copy owned by Ren Potts: both Potts and Stephenson were students of 
Cyril Domb.[H [5] It's possible that Potts's copy had come from Domb when they were 
together in Oxford from 1949 to 1951, but more likely that he had been given it by Elliott 
Montroll, with whom Potts had collaborated in the early 1960's (see ref. [6]). A transcript 
of the author's copy is given in section 3, and a scanned copy forms Appendix B. 

Both these methods start from the formulae for the pair correlation of two spins 
deep within an infinite lattice. There is also a third method used by the author for 
the superintegrable chiral Potts model (which is an iV-state generalization of the Ising 
model) [7]: if one calculates the single-spin expectation value (<ti,m) in a lattice of width L 
and height 2M with cylindrical boundary conditions and fixed-spin boundary conditions 
on the top and bottom rows, then one can write the result as a determinant of dimension 
proportional to L. (This method is similar to that of Yang. [8 1) The determinant is not 
Toeplitz, but in the limit M — > oo it is a product of Cauchy determinants, so can be 
evaluated directly for finite L.[H [TO] 



2 The first method 



In August 1948, Onsager silenced a conference at Cornell by writing on the blackboard 
an exact formula for Mn. |lH p. 457] The following year, in May 1949 at a conference 
of the International Union of Physics on statistical mechanics in Florence, Italy, [121 P- 
261] Onsager referred to the magnetization of the Ising model and announced that "B. 
Kaufman and I have recently solved" this problem. He gave the result as 



Ma 



[l-k 



2\1/ 



(2.1) 



where 



k = l/(sinh2# sinh2#') (2.2) 

and the result is true for < k < 1, when T < T c . For k > 1 the magnetization vanishes, 
i.e. A4 = 0. Figure [2] shows the resulting graph of Mq for the isotropic case H' = H. 

Onsager and Kaufman did not publish their derivation, which has led to speculation 
as to why they did not do so. The first published derivation was not until 1952, by C. N. 
Yang,[8j who later described the calculation as "the longest in my career. Full of local, 
tactical tricks, the calculation proceeded by twists and turns." [13, p. 11] 

Onsager did outline what happened in an article published in 1971. [H] He starts by 
remarking that correlations along a diagonal are particularly simple, and gives the formula 



(<7l,lCm,m) — An-1 

where D m , like A m , is an m by m determinant: 



Co ci c 2 • • • c m _i 

C-l Co Cl • • • Cm-2 
Cl-m Co 

the c r being the coefficients in the Fourier expansion 

00 



(2.3) 



(2.4) 



E 



c r e 



(2.5) 



4 



M 
1 



o Tc r 



Figure 2: M. function of temperature T . 



of the function^ 



l-ke 1L 
1-ke- 



1/2 



(2.6) 



It is clear that Onsager knew this when his paper with Kaufman was submitted in 
May 1949, because footnote 7 of [3] states that "It can be shown that 5' = ir/2 — w/2 at 
the critical temperature for correlations along a 45° diagonal of the lattice." Indeed, this 
result does fo llow immediately from (2.6) when k = 1, provided we replace 5 in (1.4) by 
5. Certainly (2.3) is true, being the special case J3 = V3 = of equations (2.4), (5.13), 



(6.10), (6.12) of Stephenson's pfaffian calculation [15] of the diagonal correlations of the 
triangular lattice Ising model. The formula (2.3) agrees with (1.7) when m = 2. 



The fact that the formula depends on H, H' only via A; is a consequence of the property 
that the diagonal transfer matrices of two models, with different values of H and H' , but 
the same value of k, commute. [161 §7-5] 

In [14J, Onsager says that he first evaluated D m in the limit m — > 00 by using gen- 
erating functions to calculate the characteristic numbers (eigenvalues) of the matrix D m 
and that this leads to an integral equation with a kernel of the form 



K(u,v) = K\(u + v) + Ki{u — v) 



(2.7) 



He then obtained the determinant by taking the product of the eigenvalues and says that 
"This was the basis for the first announcements of the result.". We show how this can be 



2 I write the 5 of HH as 8. 
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done in Appendix A. One does indeed find a kernel of the form (2.7), and go on to obtain 



lim D m = (1-fc 2 ) 1 / 4 , (2. 



m— >oo 



which is the result (2.1 ) that Onsager announced in Cornell in 1948 and Florence in 1949. 



3 The second method 

In his 1971 article [2] Onsager goes on to say that after evaluating the particular deter- 
minant -Dqo by the integral equation method, he looked for a method for the evaluation of 



a general infinite-dimensional Toeplitz determinant (2.4), with arbitrary entries c r . (The 



c r must tend to zero as r — > ±00 sufficiently fast for the sum in (2.5) to be uniformly 



convergent when u) is real.) As soon as he tried rational functions of the form 

nYl — n p iuJ \ 

m = Offl • M-lfcK 1 - (3 ' 1) 

"the general result stared me in the face. Only, before I knew what sort of conditions to 
impose on the generating function, we talked to Kakutani and Kakutani talked to Szego, 
and the mathematicians got there first." 

In another article in the same book [T7], Onsager gives further explanation of that 
comment, saying that he had found "a general formula for the evaluation of Toeplitz 
matrices]^] The only thing I did not know was how to fill out the holes in the mathematics 
and show the epsilons and the deltas and all of that". Onsager adds that six years later 
the mathematician Hirschman told him that he could readily have completed his proof 
by using a theorem of Wiener's. 

There is contemporary evidence to support these statements in the form of correspon- 
dence in 1950 between Onsager and Kaufman. There is also the photocopy of a typescript 



mentioned in the Introduction, which deals with the problem of calculating the of (1.3). 
Here I present a transcript of it, with approximately the original layout and pagination. A 
scanned copy is in Appendix B. It seems highly likely that this is Kaufman's initial draft 
of paper IV in their series of papers on Crystal Statistics (see points 3 and 4 in section 5). 

Hand-written additions are shown in red (for contemporary additions) and in blue (for 
probably later additions). Not all additions are shown. The "Fig. 1" mentioned on page 
2, after equation III. 45, may be Fig. 4 of [3J. Equation (17) on page 7 follows from eqn. 
89b of p] after interchanging H' with H* , 5' with 5* and setting H' = H. There is an 
error in the equation between (20) and (21): the second e~ 2H in the numerator should be 
e 2H , as should the first e~ 2H in the denominator. 



3 Refs. [T3], [IZ] are reprinted in Onsager's collected works, pages 232 - 241 and 37 - 45, respectively.! 
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The unpublished draft 



R.B. Potts 



B. KAUFMAN 

Long-Range order 

Onsager 

Introduction 

In a previous paper 1 we investigated the problem of short-range 
order in binary crystals of two dimensions. Short-range order in the 
neighbourhood of a fixed crystal-site (a, b) was described by a family of 
curves (s a ,&Si,j)Av. , each curve corresponding to a different choice 
of i and j. These functions give the correlation between the two crystal 
sites (a, 6) and as a function of temperature. In other words, (s 0) & s y)Av. 
is the probability that the spin at is +1 if it is known that the spin 
at (a, b) is +1. 

It was shown in III that for low temperatures 

{s a ,bSi,j)Av. ~ 1 

The physical meaning of this is clear: the model we used was a ferromagnet, 
in which, at zero temperature, all spins are aligned in the same direction; 
therefore all correlations (for all choices of ) are equal to +1. 
The correlation functions decrease as the temperature rises, first slowly, 
and then quite rapidly in the neighborhood of the critical temperature. For 
high temperatures the correlation functions tend to zero, as the crystal be- 
comes more and more disordered. (Fig?) 

Long-range order is the limit of the set of functions (s a ,& s i,i)Av. 
(for a fixed temperature) as the distance between (a, b) and tends to 
infinity. We will evaluate this limiting function for correlations within 
a row, that is to say, for the case i = a. Since only the relative positions 
of the sites matter here, we may take a = b = 1 and evaluate 

lim (si 5 iSij) A v. 

J-tOO 

It was shown in III that 
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Sij)av. = cosh 2 H*.Aj — sinh 2 H*.A_j 
Here H* is a function of the temperature given through the relations 

e~ 2H = tanh H* , H = J/kT 

J is the energy parameter of the crystal, and k is Boltzmann's constant. 

Aj and A_j are j-rowed determinants, expressed in terms of certain functions 

of temperature, S r : 



A, 



(111.43) 



(1-7,14) 



Si S2 


S3 • 


• Sj 




S_i S_2 


S_ 3 • 
















So Si 


s 2 • 






So S_i 


S_ 2 • 


• S-j+i 


S_i So 


Si • 


' S j-2 




Si So 


S_i • 


• S-j+2 


S_j +2 - 




• Si 




Sj-2 • 




• S-i 



(111.45) 



(Fig. 1) 



The temperature dependence of the S r is shown in Fig.l. 
Qualitatively we can see from the Figure that, since, for both high and low 
temperatures S_ r is very small, A_j must be small; but Aj, which has Si 
along its main diagonal, is small only at high T, whereas for low T its value 
is close to +1. 

To evaluate Aj exactly, we must know the analytic expressions for S r , 
which have been introduced in III. It was shown there that the S r are Fourier 
coefficients of a known function, exp[i5'(uj)] , and we will here make essential 
use of this fact. 



Evaluation of Aj : 



What is the value of 



A,- 



CO 


Cl 


C2 • 




C-l 


CO 


Cl • 




C-2 


C-l 


co • 




C-j+1 






• c 



(1) 



where the c n are the Fourier coefficients of a given function 
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/M = E Cn<3 

— OO 

At first we consider the special case where 



a + aie luJ + a 2 e 2luJ + ■■■ + a p e piu; 



b + bie~ luJ + b 2 e~ 2luJ H + ^e"^ h(e~ luJ ) 

and h(z) ^ for \z\ < 1 



It is no restriction to assume: 



&o = l • 



Multiplying both sides by h{e lu ) we get 



c n e 



\r=0 / V-oo / 

OO 

= E ( b o°n + hcn+i H b q c n+q ) e l 

n— — oo 

= E«* e ^ = 0(0 



t=0 



Hence 



6 c n + hcn +1 H h 6qC 







(2) 



(3) 



(4) 



for all — oo < n < oo, except the finite set: n = 0, 1, 2, . . . ,p. In particular 

b c + ha H h = a . 

(4) means, when applied to the determinant Aj (j larger than q), that linear 
combinations of rows in Aj can be chosen so that most elements in those rows 
vanish. 

For example, suppose that boc n + b\c n+ \ + 62Cn+2 = for n / 0, 1,2, 3, 

Then, in A4 = cq c\ c 2 C3 
c-i c ci c 2 
C-2 c_i c Ci 
C-3 C_2 C_i C 

multiply the 3rd row by bo, and add to it the 2nd multiplied by b\ and the first multiplied 
by 62- Similarly, multiply the forth row by bo and add to it the third 
multiplied by b\ and the second multiplied by b 2 . These operations do not affect 
the value of A3, and so we have: 
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C Cl 
C_i c 



A 4 = c ci c 2 c 3 

C_l C Cl c 2 

ao ai 
a 

All elements to the left of ao vanish, and the same will be also true for 
larger determinants. Therefore, in this case: 



A, 



,i-2 



•A, 



In general case, when h(e w ) is a polynomial of degree q, 
we will find: 



A, 



A fi 



A,- vanishes if ao = 0, that is, if there is no constant term in the polynomial 
g{e luJ ). In other words Aj = 0, if 

g (e iw ) = e iku ■ gi(e iul ) Provided j > q. If j = q can't say this! 

where gi(e luJ ) is a polynomial in e 4a; , with a non- vanishing constant term. 
Furthermore, if ao > 1, (5) shows that Aoo = lim^oo Aj diverges; 

while if ao < 1 A ra = limj^oo Aj = 0. We disregard both of these uninteresting 
cases, and consider only the case ao = 1. (As mentioned above, bo can always 
be chosen to equal 1). We can now replace (5) by 



A; 



A fi 



V > ?) 



A q can always be expressed in terms of the constants in g(e lul ) 
and h(e~ lul ). Rather than use the coefficients at,b r , we factorise the 
polynomials 



h(e~ iu) ) 



t=i 
q 



n(i-A 



r=l 



and we express A q in terms of the a's and /?'s. 
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We will now show that (1 — atPr) is a factor of A q . Consider f(co) = g{e lu) )/ 
h(e~ lu) ) as a function of the variable a r . When a r = (A) -1 we get 



g(j") 



[1 - aie iw ) •••(!- a r e 



h(e~ iuJ ) (1 - (3^^) ■■■{!- (3 t e- iuJ ) ■ ■ 

a r e iul (f3 t e~ luJ - 1)(1 - a^) 



(l-/3 i e-^)(l-/3ie-^) 



e iuJ ffi(e^) 



where gi(e luJ ), hi(e~ tuJ ) are both polynomials with constant terms. Hence 

by the above argument Aj = when a r = (Pt)^ 1 , and we conclude that 1 — a r f3t 

must be a factor of A q . Since this argument applies to all pairs a r ,(3t, 

we find 

A, = F(fii, 02,..., P q ) ■ ]~J(1 - a r (3 t ) S^sultant 



Now to find F((3i,@2, ■ ■ , Pq)- Since F is independent of the a's, we 



may choose ol\ = a2 



Otr, 



0, and then get g(e lul ) = 1 



For f{oj) we have in this case 



g(j< 



h(e~ iul ) 



assumed h{z) ^ for \z\ 



but the right-hand side will contain only terms with e -maJ ; therefore, in 
f(to) all c n = for n > 0. The determinants corresponding to this f(to) will 
have the form 



A„ = 1 •■ 
x 1 •• 
y x 1 •• 

On the other hand, our choice of a r gives 

3,p 



(x,y,...^0) 



H(l-a r pt) = 1 



As a result 



F(p 1 ,p 2 ,...,P q ) 
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and this holds for all choices of a r . 

We thus have, for a general choice of a r ,(3t- 



A q = Y[(l-a r (3 t 

r,t 

Together with (6), this gives 



q,p 



Aoo = A, = JJ(1 - a r /3 t ) 

r,t 

for the case where f(u) is the ratio of the polynomials given in (7). 
(13) can be generalized. Consider the case where several roots of 
g(e luJ ) coincide, and similarly for h(e~ lw ). Then we may write: 



and from (7) we get in this case 



nij and n k are integers, from the construction. 
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Log-range order along a row. 



In order to apply the theory developed above to the problem of 
long-range order, we have to know the function in which the S r are 
Fourier coefficients. From III we have 

1 f w 

E r = - / cos[ru> + S'(w)]dw (16) 
71 Jo 

About 5'(w) we know (from I. ) that 

. sinh2iif s'mw 

tan o (w) = — - „ rr TT 

coth 2ii — cosh 2H ■ cos u; 

<$'(u;) is therefore an odd function of w, so that we can write 

E r = 1 r e lS '^ -e trw dw 
K Jo 

and, conversely, 

oo 

e i<5'0) _ ^ ^ Jrw 



ri V 

— oo 



e tS '( w ) now plays the role of f(w) from ( ), and we have to express it as 
a ratio of two polynomials, as in ( ). For this purpose we write 



l 



e l5 ' 



1 + i tan 5' \ 2 / coth 2ff — cosh 2H • cos to + i sinh 2i? • sin w 
1 — i tan <5' / \ coth 2H — cosh 2H ■ cos u> — i sinh 2i? • sin w 



/coth2if - i e - 2ff+ ™ - ie" 2i/ -™\ 5 
I coth2# - \e~ 2H+lw - \e- 2H - lw ) 



,1-2 coth 2H ■ e~ 2H+lw + e - 4i/ + 2 ™ \ 2 

= e 



^1-2 coth 2H ■ e- 2H ~ iW + e" 4 ^" 2 ™ 
This is now almost in the desired form (3). If we define 

5(w) = 5'(w) + w 

we have 

l 

iS{w) ( 1 - 2 coth 2i7 • e" 2i/ +™ + e - 4 ^+ 2i ™ \ 2 
6 ~ V 1 - 2 c °th 2H ■ e - 2H ~ lw + e" 4 ^- 2 ™ J 
We now factor the polynomials in the denominator and numerator, and find 

l 

iS(w) f 1 - coth H . e- 2H+lw ){l - tanh H ■ e~ 2H+lw ) 1 2 
\ 1 - cothtf . e" 2 ^-™)(l - tanh H ■ e - 2 ^-™) J 
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We can now find the determinants corresponding to e lS ( w \ But we observe that 

oo 

e i8(w) _ e iS'(w) _ Q iw _ (Yj 



?? 



E 



-i)e* 



E 



Hence, 



the theory of section 2. gives us the values of 



A, 



C Ci c 2 
C-l c Ci 
C-2 C_l C 



51 So S_i • 

52 Si So • 

53 S2 Si • 



Si S2 S3 
So Si S2 
S-i Sq Si 



Aoo is the limit of Aj in (111.45). We find from (13) and (7): 10 

Aoo = {(l-coth 2 F-e- 4H )(l-tanh 2 F-e- 4H )(l-e- 4H ) 2 }2 
But from ( ) we derive tanhfP = e~ 2H , so that 

Aoo = {[l-e 4 ^- 4 ^][l-e- 4 ^*+^)(l-e- 4ff ) 2 }^ 

On the other hand, the limit of Aj in ( ) is = 0. This is so, because the 
coefficients appear in A_j belong to the function e %w ■ e l5 ^ w \ and we have seen 
that the corresponding determinant vanishes in the limit ( ). 
As a result, long-range order along a row is given by: 

lim( SllSl)i ) Av = cos^^fl-e^-^itl-e^+^Itl-e-^S 



14 



4 Summary of the draft 



The names at the top of the first page have been added by hand to the typescript, so it 
is not immediately obvious who are the authors. However, from the first sentence and 
their frequent references to paper III , specifically to III. 45 (particularly in conjunction 
with their correspondence discussed below) it is clear that the paper is by Onsager and 
Kaufman jointly. Footnote 1 in the first line is not given, but is presumably their paper 
III. 

The paper begins by quoting III. 43, modified to (1.2) above, for the isotropic case 
H' = H. 

It focusses on the problem of calculating a j by j Toeplitz determinant A,-, of the 
general form (2.4), in the limit j —> oo. The c, are the coefficients of the Fourier expansion 
(2.5) of some function f(oo), initially allowed to be arbitrary. It first takes f{oo) 

and shows that 



q + aie^ + . . . + a p eP iuj 
b + &ie~ iw + . . . + b q e-i iui 



A,- = if ao = and j > 



It then takes f(oj) to be of the form (3.1), or more specifically 

5 (e-) 



where 



g(j») = H(l - a t e k 
t=l 



/i(e" iw ) ' 



(4.1) 
(4.2) 

(4.3) 



r=l 



and goes on to show, using (4.2), that 



p i 



i=lr=l 



provided that j > q. This is an algebraic identity, true for all at,/3 r . It is a trivial 
generalization then to say that if 



then 



A, 



n^ = i(i-^ 
nLia-Ae- 1 



p q 



(4.4) 



j=l k=l 

for positive integers ra^n^. 

For Aj, f{oj) = e lS ^ is given by (1.6). This is of the general form (4.4), withp = q = 2 
and 

di = /?l = coth-ffe -2 ^ , 02 = 02 = tanhi/e -2 ^ , mi = mi = n\= n<i = \ 
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so the mj,rij are no longer integers. The paper assumes that (4.4) can be generalized to 



this easel 4 ! so obtains (using tanh H = e~ 2H * ) 



{[1 



■*)][! 



-A(H* 



+H) ]{1 



-4__j2jl/4 



(4.5) 



If we transpose the second determinant ( |1.3[ ) in ( |1.2[), then its generating function is 
not /(w) but e 2lU3 f(uj). This corresponds to the form (|4.1[), but with oq = (and a\ = 0), 



so A-oo should vanish because of (4.2). Thus (1.2) gives 



lim (anan+j) 

j-.oo 



cosh 2 H*{[1- e 4 ^*-")] [1 - e- 4( - H * + V][l - e" 4 "] 2 } 1 / 4 . 



For the isotropic case H' = H, e 2H = tanhiif, so we obtain, using (1.8), 
r2 - lim (<7i,i(7i )1+j ) = (1- l/sinh 4 2^) 1 / 4 = (1 - k 2 ) 1 ^ 



Aft 



J-KX) 



in agreement with (2.1). 



Generalization to the anisotropic case 

In their paper III, [3j Kaufman and Onsager focus on the isotropic case H' = H, as does 
the above draft. It does in fact appear that their results generalize immediately to the 
anisotropic case when H, H' are independent, and H* , 5' are defined as in [Tj (i.e. by eqns. 
(1.1), (1.5), (1.6) above), and that (1.2), (1.7) are then correct as written^] Replacing 
equation (24) of the draft by (1.6), and repeating the last few steps, one obtains 



lim (a^a^x+j) = cosh 2 H* {[1 - coth 2 H' e - iH ][l - tanh 2 H'e~ 4H ][l - e" 4H ] 2 } 1/4 



[1 - l/(sinh2i/sinh2 J H" / ) 2 ] 1/4 = (1 - k 2 ) 1/4 , 



(4.6) 



which again agrees with (2.1). 



5 Further comments 



There are several letters on the Onsager archive in Norway between Onsager and Kaufman 
relating to this second method, at 



http: // www.ntnu.no /ub /spesialsamlingene / tekark / tek5 / research/009 _0097.html 



and 



http://www.ntnu.no/ub/spesialsamlingene/tckark/tek5/rcsearch/009_0096.html 



We shall refer to these two files as 0097 and 0096, respectively. 

1. On pages 21 - 24 of 0097, (also in [19j ) is a letter dated April 12, 1950 from 
Onsager to Kaufman giving the argument of the above draft up to equation 15 therein. 
Onsager states that the result admits considerable generalization. He defines 



r ]+ (oj) = logg(e i n , t?_ (w) = - log h(e~^) 



(5.1) 



4 This vital point is considered by Onsager in the first letter quoted in section 5. 

5 I have verified that (1.7) agrees with (2.3) - (2.6). I have not verified (1.2) directly, but have checked 



algebraically, with the aid of Mathematica, that it agrees with eqn. 56 of [5] for j < 
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so /(w) = e v+ ^' )+ri ^^ UJ \ He remarks that if log f(uj) is analytic in a strip which contains the 
real axis, then these functions may be approximated by polynomials which have "no wrong 
zeros" in such a manner that the corresponding determinants converge. This implies that 
the otj, (3j all have modulus less than one. He goes on to say that (15) is equivalent to 

i f 27T 

logAoo = — / r ?+ dr ? _(a;) , (5.2) 

If 

oo 



and bo = 0, then (5.2) is equivalent to 

oo 

logAoo = ^nbnb-n , (5.4) 
n=l 

which is the result now known as Szego's theorem. [20 |, 1211 [6] 

Onsager then says "we get the degree of order from C.S. Ill without much trouble. It 
equals (1 — k 2 ) 1 ^ 8 as before." 

2. On pages 32, 33 of 0096 is a letter from Kaufman to Onsager. It is undated, but 
was presumably written after April 12th and before April 18th 1950 (the date of the next 
letter we discuss). She thanks him for his letter and comments that his method is elegant 
and simple, far superior to the integral-equations method for this purpose. 



She also mentions the formula (2.6), i.e. 

1/2 



1 - ke 11 
1 - jfce- 



for the generating function for long-range order along a diagonal. 

She goes on to say that the mathematician Kakutani had written to her saying that 
he had spoken to Onsager about this. He had been very interested in Onsager's letter, 
which he saw in her house and immediately copied down. 

3. On pages 30, 31 of 0096 is a letter from Kaufman to Onsager dated April 18, 1950, 
she says she's glad to hear that Onsager is coming to Princeton the following Friday and 
refers to Onsager's 1944 paper I, to her 1949 paper II, to their joint paper III, and then, 
significantly to a paper IV "on long-range order" . She says she would like Onsager to see 
her m.s. 

4. Then on page 34 of 0096 is a letter from Kaufman to Onsager dated May 10, 
saying "Here is a draft of Crystal Statistics IV" . 

All this fits with Onsager's recollections of 1971 |14l I17j. He and Kaufman had evalu- 



ated the determinants in III.43[3], or alternatively in (2.4) - (2.6), by the integral equation 



method (Appendix A) and by the general Toeplitz determinant formulae (15) of the draft, 



i.e (5.2). He preferred the second method, but apparently was content to let Kakutani 
and Szego take it over and put in the rigorous mathematics. Kaufman's draft of paper IV 
was never published, but it seems highly likely that the draft given here is indeed that. 

Szego did publish his resulting general theorem |20j. |2H p. 76] on the large-size limit of 
a Toeplitz determinant, but not until 1952. He also restricted his attention to Hermitian 
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forms, where f*(w) = f(oj) and b* n = 6_ n ,P, footnote 17] so his result needed further 
generalization before it could be applied to to the Ising model magnetization. 

The first derivation of Mq published was in 1952 by C.N. Yang.|8j He used the spinor 
operator algebra to write Mq as the determinant of an L by L matrix and evaluated the 
determinant by calculating the eigenvalues of the matrix in the limit L — > oo. Intriguingly, 
he mentions Onsager and Kaufman's papers I, II and III, and in footnote 10 of his paper, 
Yang thanks Kaufman for showing him her notes on Onsager's work. However, his method 
is quite different from theirs. 

Later, combinatorial ways were found of writing the partition function of the Ising 
model on a finite lattice directly as a determinant or a pfaffian (the square root of an 
anti-symmetric determinant) . [22j 123] . Then it was realised that the problem could be 
solved by first expressing it as one of filling a planar lattice with dimers. [24" ] 125 ] 126 ] 127] In 
1963 Montroll, Potts and Ward[6] used these pfaffian methods to show that (cr^icij+i) 
could be written as a single Toeplitz determinant. They evaluated its value in the limit 
j — > oo limit by using Szego's theorem. 

So there is no reason to doubt that Onsager and Kaufman had derived the formula 
(2.1 ) by May 1949, and ample evidence that they had obtained the result by what is now 
known as Szego's theorem by May 1950. They did not publish the calculation, perhaps 
because the mathematicians beat them to the remaining problem of "how to fill out the 
holes in the mathematics and show the epsilons and the deltas and all of that" . 
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Appendix A: An integral equation method 



We regard D = D m as a matrix and write z for e lw and C(z) for e lS ( u \ so (2.5), (2.6) 
become 

1-kzV' 2 ^ 



c « = [T=kT> = 1* c "' ■ <A1) 

N ' / m=—oo 

For < k < 1, the determinant of D m tends to a non-zero limit as m — > oo. The 
eigenvalues of D itself lie on an arc in the complex plane, and appear to tend to a 
continuous distribution as m — > oo. By contrast, in this limit the individual eigenvalues 
X r of D T D occur in discrete pairs, lying on the real axis, between and 1. If one orders 
them so that A r < A r +i, then for given r the eigenvalue A r appears to tend to a limit as 
m — > oo, and these limiting values tend to one as r — > oo. 

We assume these properties and seek to calculate the eigenvalues of D D in the limit 
m — > oo, and hence the determinant of D T D, which is (det-D) 2 . Writing them as A 2 , we 
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can write the eigenvalue equation as 



\x = Dy , Xy = D T x . (A2) 

For finite m, x = {x Q ,x 1 , . . . ,x m -i}, y = {yo,Vi, ■ ■ - ,2/m-l}- Let P be the m by m 
matrix with entries 

P%j — 3i,m—l—j 

so Px = {x m -i,x m - 2 , ...,xq}. Then 

D T = PDP , D T D = (PD) 2 . 

For finite m, the eigenvalues of D T D are distinct, so the eigenvectors are those of PD. 
For large m and a given A, this means that the elements Xi,yi are of order one if i is close 
to zero or to m — 1, but tend to zero in between, when both i and m — i become large. 

However, in the limit m — > oo, the eigenvalues occur in equal pairs. One can then 
choose the two corresponding eigenvectors so that one has the property that 

Xi,yi^-0 as i — > oo , (A3) 

while for the other eigenvector x m _j, y m -% — > 0. The eigenvectors are transformed one to 
another by replacing Xi and by y m _i_j and x m _i_j, respectively. 

Since the eigenvalues are equal, we can and do restrict our attention to eigenvectors 



with the property (A3). 



Generating functions 



Taking the limit m — > oo, we can write (A2) explicitly as 



Xxi = '^2c j - i y j , Xyi = '^2c i - j Xj (A4) 

3=0 3=0 

for i = 0, 1, 2, ... . 

We can extend these equations to negative i, defining Xi,yi to then be given by the 
left-hand sides of the equations. Let X(z),Y{z) be the generating functions 

oo oo 

X(z) = ]T xiz 1 , Y(z) = Hi** ■ 

i=— oo i=— oo 



Then (A4) gives 

XX(z) = C(l/z)Y(z) , XY{z) = C(z)X(z) , (A5) 

where 

oo oo 



X(z) = Y J *iZ i , Y(z) = J2vi 



A Wiener~Hopf argument (or simple contour integration on each term) gives, for 



X(z) = <b dw , Y(z) = <p dw , (A6) 

w 2vri J w - z ' w 2vri J w - z ' y J 

the integrations being round the unit circle in the complex w-plane. 



Thus (A5) is a pair of coupled integral equations for X, X(z),Y(z). 
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The integral equations in terms of elliptic functions 



To get rid of the square root in (Al) we introduce Jacobi elliptic functions of modulus k 
and set 

z = k sn 2 u . 

If K,K' are the complete elliptic integrals and u = a — K — iK'/2, then from (8.181), 
(8.191) of [28], or (15.1.5), (15.1.6) of HB|, 

~ (l +p 4n-3/ s )2( 1+:p 4n-l s )2 

Z ~ S 11 (l+p4n-3 5 )2( 1+p 4n-l/ s )2 
n=l 

where 

p = 9 l/2 = e -nK'/2K ? s = j*«/K j < g < 1 . 

If < k < 1, then < p, q < 1, so we see that z goes round the unit circle as a goes 
from —K to if along the real axis, i.e. u goes along a horizontal line in the complex plane 
from -2K - \K' /2 to -iK'/2. 

Prom the elliptic function relations cn 2 n = 1 — sn 2 u, dn 2 u = 1 — k 2 sn 2 u, 

„, . . snndnii 

C(z) = i , 

cnu 

tV (l + (-lTp 2 ^s)(l-(-lYp 2 ^/ S ) 

11 (i-(-i)V n ~ 1 s)(i + (-i)V n - 1 A) 



Replace w in (A6) by 

w = k sn 2 v , 

where 

Im(v) = -if/2 , -iT'/2 < Im(u) < K'/2 . (A8) 

Then 

dw = 2k sn v cn u dn v dv . (A9) 

It is helpful to set 

/ \ —iX(u) , . iy(tt) 

X « = -p- , Y(u) = . A10 

sn u dn u cn u 



Then the integral equations ( A5 ) become 

\X(u) = — [ M(u,v)Y(v)dv , XY(u) = — [ M(v,u)X(v)dv , (All) 
2vr J 2ir J 

where the integrations are from — 2iT — IK' /2 to —\K'/2 and 

, r , . „ cn u sn v dn v 
M(u,v) -- 



sn 2 v — sn 2 u 



Using the Liouville theorem-type arguments of sections 15.3, 15.4 of [16], one can establish 
that 

. ,/ x dn(v + u) dn(v-u) 

M(u,v) = — ) r + — 7 f • 

sn(u + n) sn(v — u) 

We see that this kernel is indeed of the form (|2.7 
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We require the functions X(z), Y(z ) t o be anal ytic fo r \z\ < 1, i.e. for |Im(u)| < K'/2, 
and in particular for real u. From (A5), (A7) and (A10), X(z) = XY(u)/snu dnu, Y{z) = 
XX (u) /cnu. Since snO = cnK = 0, it follows that 



X{K) = Y(0) = . 



(A12) 



Also, M(u, v) is an analytic function of u within the domain ( A8 ), and we can negate u 



within the domain. Since M(u, v) is an even function of u, it follows from (All ) (provided 



A is not zero) that X(u) is also an even function. Similarly, Y(u) is an odd function. 



Solution by Fourier transforms 



The sum and difference form of the kernel M{a,f5) suggests solving (All) by Fourier 
transforms. The functions X, Y, M are anti-periodic of period 2K, while X(u), Y(u) are 
even and odd, respectively. We therefore try 



X{u) = ^cosvr(2r- l)u/2K , Y{u) = £>sin7r(2r — l)u/2K , 



(A13) 



r being a positive integer. We note that this ansatz immediately satisfies (A12). 



Integrating round a period rectangle and using Cauchy's theorem, or using (8.146) of 
, for all integers r we obtain 



1 

2^ 



dll (^ ~ U ) ^(2r-l)v/2K dv 
sn(v — u) 



j-n(2r-l)u/2K 



1 + 



■,2r-l 



(A14) 



Negating 2r — 1 and/or u (but not v ) and taking appropriate sums and differences, we 
obtain the identities, true when both vq + u and vq — u have imaginary parts between 
—2K' and zero and the integration is from vq — 2K to vq : 



1 

27 
1 

2vr 



M(u, v) sin 



n(2r-l)v A 
2K aU 



1 



1 + 



I n(2r-l)u 

2r _i cos 27? 



M(v, u) cos ' K ^ I ^ V dv 



1 



2r— 1 

Q ■ 7r(2r-l)tt 

1 + Sm ^7^ 



We see that the ansatz ( |A13 ) does indeed satisfy (All), with A = B = 1 and 



„2r- 1 



1 + 9 



2r-l 



(A15) 
(A16) 

(A17) 



More generally, we can allow X(u), Y{u) each to be an arbitrary linear combination of 
exp[i7r(2r — l)u/2K] and exp[— i-7r(2r — V)u/2K\. We obtain the above solution, plus an- 
other where A = 0. However, this second solution does not satisfy the necessary condition 



(A12). 



Determinant of D 

For every positive integer r there is one and only one eigenvector of D T D of the form 



(A3), with eigenvalue A 2 , A being given by (A17). However, from the discussion before 
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( A3 ) , there is also an equal eigenvalue whose eigenvector has elements re, , yi of order 1 
only when i is close to m. Hence, using (8.197.4) of |28] or (15.1.4b) of |16j . 



det 



D T D = (detDf = H(^ 



x / l_^2r-l\4 



r=l 



+ 



,2r-l 



k> = (1 - k 2 ) 1 / 2 



so 



det D = (1-k 2 ) 1 ^ 



(A18) 
(A19) 



as in (2.8). 
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Appendix B: The typewritten draft 



ft 8 .. Po^ 



Introduction 



In a previous pcpsr^ ttb Irrrestlcatod Ve problem of short-* rftuvg^ 

order in binary crystals of fcpo distensions, S^ort-ren^e oilier In .the 

• nelgbboxiiood of a fixed crystal-aite (£>b) ttss described by a fa-ally c£ 

curares i^" tr, , each co-ye corrpsjxjadins to a different c-&c-ic~e 

of 1 imd ^ * Ttese functions give tlx correletion between the two c^jibal 

oitoa (a,b) and (1, j) as a function of te^xmyturo. In other -wordo J <'a a J 

is the probability that t:ie .spin at is ->1 if it is knom that taie opin 

at (a,b) is +1* 

It was sfiorsn in III tnat for low t^per&tnres 

<fa .s, />. 
v a, b i, j / Av. 

The physical aoaning of tide is clear i the ssodsl we esed -san a ferroosgaet, 
in wbdefaj at sero teraperature, nil spine -are aligned in tha satae dirBCtioaj 
therefore. All co reflations (for all choices of (i>J) ) &re equal to 
The correlation functions decrsase as the tewperatore rises,, first aLoTsly, 
and then -quitfl rapidly in the neighborhood of t?,© critical t&a'per&tare* For 
high temperatures the correlation functions tcx:d to sero, ae the crystal be— 
costsa nor© and isore disoi-dtered. (Tig?) 

Locxf-itsTft-a- ordsr is" the liait of t'a pat of functions*^ 8 , S , .N„„ 
— iss ~& — — ^ a # b i^J ' AY. 

(for a fixed temperature) as tl»e dietaries botsrsfin (a,b) and (I* J) tends to 
infinity, ";e -will ovslmvte this Halting function for correlations -Kithin 
a" rem, that is to gny, for the oa.se i — a, Since only the relative positions 
of tha' sites astter here* "we nay take a ** b c 1 end evaluate 

It ^as skotia in IH tiist 

■.'.V -.• 23 " , 



"G : >> ; » > t ^ t - coc'tSw. \ j - sixth?. A.~S (III, Ls3) 

A " - i <i / J|. - * * * 

<} 

Ibi.? L 1p a fraction of the tc^psratarci given through tha relations 

eT 2 ^ toil / » H = J/kT (I. -7 , 14) 

J is 'uho enarjy par&rratcr of the crystal, and k is Boltisaann' o constant. 
/\ ^ and . ■^-„^ ar * j- rowed det3rrainsntB, oppressed in terras of certain functlorA. T 
o£ teapsrature, Z„ s 

A-f 



1 H 3 j 



— 1 O i. 



3-2 



f> C « * A 



Z 



-1 <-t -j -J 



■j * 1 



i 

1 * * « * » 

Si, 



>2' 



(Fig, D 



Tisa temperature dependence of the 2^ ia shomi In Fig. 1. 
Qualitatively we can see from tha Figure that, since, for both high and low 
tsaperafcuras Z r is vsry sasll, ^ _^ aust be saallj but A^, which baa 2^ 
along it© main diagonal, is saall only at high T, whereas for low T its value 
is close to + 1. 

To evaluate exactly, yie must know the analytic expressions for 

which have been introduced in III. It was shown there that the 2 are Fourier 

r 

coefficients of a known function, exp [i £( , and we will here raaka easentiel 

ues of this fact. 
Evalu ation of f 
trhat is tha value of 



c , c 
-1 o 



c -2 c -l 



-j + 1 



e a * « * 



(i) 



whers the Fourier coeffienfcs of a given function 
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» 3 - 



At first Ti\i consider the spc-sial c&so wh' ro 



(2) 



Hi?) 



b b,o 
o 1 



It la Qj> restriction to aesuaai 



b « 1 
o 



-] 



-iw, 



tisliiipJylRs loth sides by h(e ) we Get 



~c-3 n 



f (bo + b,c ..4 . . * + b c, ) 
o n 1 n+1 q it+Q 



in" 



h -- ■ 
P 



- u a t e * e (a ) 



Hence be + fc,c .-, + .*» + bc, =0 
on In* 1 qn+q 



(3) 



for ell • 03 <^n <C 00 , exespt the finite sets n = 0,1,2,, , , ,p. In particular 



b o + ko, •*■».♦ + b c = a . 
o o 1 1 q q o 

(U) scans, xhen applied to the doterolnant A ^ (J larger than q), that linear 

combinations of rows in ^ can I*) chossn so that aost p.lasients in thoas rows 

rani eh. 

For esKraple, suppose that b^c^-*- b i c n«ri + ^2°n+2 = ^ for n ^ 



Than, in E 



°o c l °2 c 3 
c ~l°o <T c 2 



c „ o ,. c c. 

— £ — 1 o 1 



c -3 C -2 G -l C o 



mltiply the 3rd ros by b^, and add to it the 2nd arultipliad by b^ and the first 
ra\xltiplied by bg. Sinilarly, multiply the forte row by b Q sM add to it the third 
aoltiplisd by b^ snd the second multiplied by b^. These operations do not af f e-o t 



the valtia of-A^. aM so «a have j 
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~ k - 



°-X % °1 °2 

i), 

a 



2 | 
° / 



o 1 

c , c 
«1 o 



Ail elements to the left of a venieh, end the earns will be also true for 



l*rc*i" <2et.crsLinaE.t-c* Therefore, in this case: 



& » a J 



J 



In gonsral case, when h(e ) is a polynomial of dsgres q ; 



we rill flndt 



-r o 



(?) 



A vanishes if a « 0, that Is, if there is no constant tern in the pol^notsial 



_iAJ 



g(e )* In other words ^ c 0, if 



iw 



: at, <■ 1^7 



«her® SkCS. } is a polynomial ia a , with a r»n~vanifth±ng constant term* 



Furfchornoro, if a > 1, (?) shous that 

0>> ?£ 



A TO = 11m /s. divargsaj 



while if a< 1 



lira ^ , « 0. 



J' 

Ve dierec ar ^ of these ualnta resting 



cases* and consider only the case a Q c 1» (As mentiori0d above, b Q can always 



be choeaa to equal 1)» can now replace (5) by 



(6) 



cm be expressed in terms of the constant© in g£e^" A " > ) 



and h(e"^"'^), Esther than use the coefficients a^, b^j, we factoriss th© 
polymslEls 

g^ itJ ) = ^ (1 -^e 1 ") (7) 

/pi J r 

and we espresa -A ^ in terras of ths-c*s andj^s* 



26 



will wrs f.ha-w tb.&tfl j^is a factor of Coiicidsr l(O) » s(^ r )/ 

L(cT^ W ) as a function of tho variable . f:hcn^' » ( .3* we rot 

r r y T. 

il^J. ) .... (I - IfLjg ) 



• * * • 



whcrs E f (^ a> )i M*""* - " ) sra V-d'Ji pelynoBiale rith constant terse* Hence 
by tha above erguacst >A ^ ■ T.he.v-S.r- ty^)""^' 188 concluc ^ 3 that 1 
Ruct bo a fnctor of A. , .Since this a^pracnt applies to sll pnirs*^, fb^ t 
re find 'C^HvJ: 

I'o-.j to find i{J2>*j% t • • • 1 £ q )» Sinca F is independent of the <Ve, wa 
choose "< / 1,r *^= . . . E x p = 0, end then L:(ti- i<U )= 1 

For £( (X/) hcve. in f'dc cage ' 

v ' . htc--^) . [ * > ITi' 1 " 

but the rij -irt- h.m<l rice wi ll _contc' -in only ter.is r.lth e"^ n therefore, in 
f(£«j ) ell c » fcr n>0. Tno £et«minrnto oorreeponlLr^ to this will 



hr.vc- the for.a 



A =,1 . . . . !" 1 

' A «L v * * • 



., . :. ■■. .:: :.: 'As: 



(::, 7, . . . *' 0) 



(IC) 



f • ■ : ... H i; 
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t 

t: 1 fj.c> bulds aw for sOl ehoicco of<^ , 

V*s thuc hevt, for a gsneral choice of ^ ^ Jh± 1 
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Q r t t v r/ V 



TOG'sthe.T *rith (6), this £ives 

4p'A^-^ d-^ft) (13) 
tor tha saera whera ) is the ratio of ths polynosnials eiircn in (7), 
(jj"3 ) csji ba g gnarsl issd. Coaoidsr the case where several roots of 
eU 1 ^) coincide, acd adLsailmOy for hCe" 1 ^. Rich we may write t 

and fros (7^) tra get in tiis casa 
S ,J ^ °k ars integers, from the construction* 



(lit) 



j\j j - 1.1 
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Int.'; rcB&g order clcru; n rva. 

la ore'.-/.' to rvply the thtory developed cbo?o to ths problfe-a of 

lai^-res^a otovv, ws have to know" ths function in which the 2 &rv 

r 

l-buricr coefficient!:* Free III havo 

*\J" oosfw * ^'(w)]dw .'^ (HI.) 

c 



Abotrfc iT'C^) tso taso-s (froa I. ) that 

fx ffM - Si^K^lnj,^^^^ = • ^ 4 T k 

v cotb £H - eochzH.coeW " c ,K-.i'-u. 

•Crr) is therefore an odd function of -sr f so that we osn write 
. „ 1 /"* irf"t( w ) ir * . v 

antf, convsreely, 

_o» — r v t- 

g 5 " ^ now plays the role of f(w) fros ( }, and we have to express it 

a ratio of t?jo polynomials, as in ( )» For this purpose ttb writs 

J- ^ 1 m A * i tan „ , coth 2H - cosh 2Hg- c og tt + i einh 2H» Sim/,. 

^T^TTTHnTT^ ~ cotETTS - cosh & ~ cos « - , raKTTir5S' / w 

iv 5 „ 1 -2H + iw 1 - 2K - iw 1 
^ cotfi 2a — -j e — e $r 

' cothSH-l^^^-le— ^ r ^ > ' 

This is uoar al-.ost in the desired form (3 ), If we define 
f(n) » £(w) + w — - 



we have J 



s 



now factor the polynorsiala in the denominator and rranserstor, and find 
e i H") j( l ~ cote H. 3 - 2K4 ^(L-tonb H . gjf 4 

((1 - coth a ."^ ~ lw )(l - tanh H . ~ ^ 
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*- ? Ot.ii ts-v; iln-i" the dsisasinants corresponding to e" 



1^ ! 



, But observe that 



rf-«* -or 



I tile theory or section 2* £Lves ^ the Tains of 

*1 Z o 2 ~1 * * 



L d 



* ft * * 



**0 "1 * * * , 

f 2 1 o j 
^ ^3 2 2 Z l * * 



• «»*** 



2 i 2 2 z r * * 



* « ♦ * s * ♦ 



A to is the Halt- of A. in Rs find froa (#) and (-.)» 



10 



^ - j(l - cothS.e'^Xl - tanh 2 !! . ^(1 - e ^) 2 



1 



But frsss, { } derive tstnb H tt « , ao that 



On the other hand, the limit of in ( ) ie « 0. Ihia is so, because tha 
coeficiont? appear in belong to the fraction B^e^ 1 ^, end w hays e**n 
that the corresponding deteEainsnt vanishes in the limit* ( }» 
is a rocult, long-rang® order along a row is giren hyi 
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